
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



THE 

AMERICAN 
MATHEMATICAL MONTHLY. 

Entered at the Post-offloe at Springfield, Missouri, as Second-class Mail Matter. 

Vol. IV. FEBRUARY, 1897. No. 2. 

EMPIRICAL FORMULAE FOR APPROXIMATE COMPUTATION. 



By the Late ANSEL N. KELLOGG, of New York City. 



[Note. The following paper is here printed in the exact form in which it was left by the author at 
the time of his death, except only a few necessary verbal alterations which are distinguished from the 
other parts of the paper by an enclosure of square brackets.] 

The following symbols have the same signification, throughout these form- 
ulae, the only exception being when the letter "d" is used with a, b, c, etc. 

n = — =index of root. 
m 

m = — =index of power. 
n 

7=quantity, or — , if fractional. 

s=sum of q-\-j> (or the semi-sum).* 

<Z= difference of q and p (or the demi-difference).* 

t- 8 . 

J7=the sub-square, or under-square=d/<)>— l) 2 . 
t/ 4 =the under-fourth=(i/</— l) 4 , etc. 
.E=Napierian logarithm of number. 
/if=]ogarithm of number to base 4. 
r=root of number. 

*These last values for s and d cannot be used in the same equation with p and q. 
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=1= signifies nearly equal. 

=||= signifies very nearly equal. 

Mercator's formula for the extraction of roots of numbers near unity was 
equivalent to my own formula No. 1, as shown in Hutton's Tracts on Mathemat- 
ics, Vol. I. 

^Z=*V q, nearly .(1). 

Hutton says he gave a formula for the correction of this result, but I have 
never been able to find it. Hutton himself gives the derivation of the above 
formula. 

First correction of above : 

,. (w* -l)rf 2 
ns + d s— - 

» ' _ * ns 

V <H (n*-V)d* 

ns—d s 

3ns 

w 2 — 1 
Substituting t for — s — , we have 

, td* 
ns+d 

V^H g (2). 

ns—d — 

MS 



Second correction of above : 



J-2- =1= WS . f gx 

\ v t(A% j.r./,._./»i-HV W # 



M8 

or by reducing, 



„ . (n»+2)8+(2w«-2 ) 1 /g»-rf»+3wd 



g=l=- — /y — (4). 

(n»+2)8+.(2n*-2)|/«*-d*-3nd 



This nearly equals 



, , w 2 -l d 2 

v<h- : . 2 -i *? < 5 >- 

ns-d—^-x-—^ 
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, 


2tU 




n 


ns—d — 


2tU 
n 



Simpler than these, and nearly related to (3) is 
ns+d ___ 

MS 

v«-- A td r (6). 

ns 

ns 

Very simple and excellent is 

ns- 
V«-<- oTfr (7). 



All the foregoing are what I call equidistant processes, because for all val- 
ues of n, the difference between the numerator and denominator of the result is 
2d, (or some multiple) ; that is, the subtractive corrections following d are the 
same in both terms, whether q/p be a proper or an improper fraction. 

Of the same nature, but differently derived, is formula (8) and its equiva- 
lent (8 J). 

■ <2«» -<*») + d»+2d» 
Vq n{2s s -d !i ) + d 2 -2ds l ;- 



. >„ , nVs* -d*)+d(2s+d) 
'' q ~ l ~n(2s^-d')-d(2s-d) mh 

The fact that all these fractional formulae are symmetrical makes their 
operation comparatively simple. 

In extracting roots of high numbers 



m n I a 



and we are thus always enabled to use q/2 m between the limits i and 2. Hence, 
the following equation becomes of value : 



„ , n'—m 
iln + m s — 

n*— m* 



nlm V^ ^Z^> ( 9 ). 



3n—m— 



9n 
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m 8 — m 2 



3'n + wi- 



9n- %i ~ mi 



3n 

or »/™,/2=l= 1 ■ , , (10), 

1 „ to*— m* 

dn—m — 



S)n 



to 2 — m 2 



3w + m— 
3n—m 



'Sn 
29(n s -m a ) 



507to 
29(ra 2 -m 8 ) 



or »/ m ,/2H= w^br W- 



507n 
In a latent form, the equidistant principle is also present in the following : 

" \l , (2w 8 +3w + l)^+(8>i 8 -2)^ + (2w a -3w+l)y 2 fl „ 

\p ' (2n 8 -3n + l)<2 8 +(8w 8 -2)cy) + (27i 8 +3n + l)p 8 l J ' 

In some of the following applications p is taken =1. 
This [viz. (13)] becomes : 

fo *v, -$±jg±£ (1 3- 2) 

f -V, JJ^g| (13-3) 

for\/-i- 15q 8 + 42qp + 7p 8 

(13-5) 



for', q 11^+339 + 6 

,? 6<j 2 + 33<j + ll 

for°,/o 91q 8 +286q + 55 



for 



\ p 13q 3 +65qp + 20 K > 



for»i/o 51<2 8 + 170q + 35 

f0r l/<2 35^+T709ToT (13 ~ 8) 

e 9 / 95q s +323q + 68 ... m 

for V<? 68^+3239 + 95 (13 ~ 9 > 
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for 10 /« 77^ + 266^ + 57^ n3_irn 

f0r \y 67q» + 266^+77iT* (13-10). 

for ««, o J767«»+a666 fr+6B67 (18-100). 

1 ,q 6567q* + 26666^+6767 V ; 

For very high indices use, without sensible error, 

" ll ... (2n + 'S)q i +»npq + (2n-S)p a - 
\p (2n-Z)q i +»npq + (,2n + '6)p ls '' 

(14) is the equivalent of (2). 

Upon the logarithmic function depend the following formulae : 



, E IE r-\\ 

„ / 2 \6 r+1/ ,,.. 

Vfl"'- g 7^ ^ZTT < 15 )- 



(E_ r—l\ 
U X r+J 



3 + ^+ ^ 

V^JL^L (16). 

3~— 

J?, the logarithm ofq-t- ( *ffijT+2 9 ( 17 )> 

3(o s — i) 8 ) 

or, if o be fractional, =|=; — , - ^ , ■ (18), 

1 (q+y) 2 + 2qp 

or, in terms of rf and 8s °l=r8~» (19)- 

This value of E, if q be between .9 and 1.1, is true to the seventh decimal, 

but may be corrected with very great accuracy, even up to the ninth or tenth dec- 

4d 8 
imal, by adding to the result the ■■ . _ a th part of itself. 



45s 



If, however q be so great as 



1.7 or so small as .6 use 44 in place of 45 

1.8 or so small as .56 use 43 in place of 45 

1.9 or so small as .53 use 42 in place of 45 
2. or so small as .5 use 41 in place of 45. 
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The accuracy of these formulae will appear from the natural logarithms 
in the next [paragraph]. 

[Some verifications of formula (18) : 
FoTq/p—U:— 

By the formula : log(^/p)=^ r =.0953101, 

error— 1 in last place. 

By the tables : 2.3978951 

2.3025851 



iog(q/p)= .0953100 
For 4/p=}f .— 

By the formula : log(q/jt.)= T l} f =.0645385, 

correct to last place. 

By the tables : 2.7725887 

2.7080502 



.0645385 
For q/p=\U •— 

By the formula : log(q/p)= r 8M, =.00995033.] 
Reconverting this by (15) we have 

1 + .004975165 + .00000825 1 _ 1.004983416 
1-.004975165+ .000008251"" .995033086' 

To this denominator add the 100th part, 

.995033086 
.009950331 



1.004983417 



Hence the real number is J#J. 

From the foregoing we have another value of "y/q as follows : 

*»-\ d -V+l){ q * + 4q p+p t ) 

(Of course r can only be taken crudely, but may by successive steps, until 
the required approximation is reached.) 

Another formula akin to (15), for values of q (or r) in terms of E is 

, E , E* 

W+ ~2~ + T2 
.,.,,_,_____ (15J ). 

n ~2- + 12 
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Let us call c s which is equal to *-=— , a root centre, meaning thereby that, 

c+2 Ti 

for all values of q, and degrees of n, — ^ nearly equals n \/q. Then 



2 + 6(r+l) 

,._1 f„_1Vr-1> ~ l_ ^ 9 U1) " 

me- 



q-1 (q-lXr-1) 

16 2~ + 6(r+l) 



2i/?+ 



9 + 1 



e, if q is under 10, is nearly ^ • (22). 

o 

c, if q is between 50 and 250, is nearly same ~~ (23). 

c, if q is between 10 and 50 [correction not given]. 

Factor Process. 

Take — x -r X . . . x — =q ; then — =q and the series is consecu- 

a b 9 « 



tive. 



„ , b d h „ bd h 

Take — x — X X — =q; then -=q. 

a c 9 «« ■ • • 9 



Now if, to take the wth root of q, we assume n terms, consecutive or non- 
consecutive, and nearly equal in value 

Wq== vi — v i VI (24)) 

,/?> ^ v 'd + + i/fc 

wherein the numerators are the square roots of the terms and the denominators 
the reciprocals thereof. 

If the terms are consecutive and odd in number 

ah + 2bg+ +2de _ 

vq 2ag + 2bf+ d* ' ^ &) 

but if consecutive and even in number 
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»,/.-■ _jgj±j ^+ ^Ai± e2 (25i). 

1/(2 2a/i + 2fe 3 + + 2de ^ * J 

Also if the series is consecutive, 

„ , = (g + h)b + (f+g)c+ +( h + c)g + (a + b)h 

' q (g + h)a + (f+g)b+ + (ft + c)/+ (a + b)h - 

(25) and (25J) are called the diagonal process. 

To separate a quantity q/p into consecutive factors which shall be nearly 
equal, and which shall be as many as there are units in the index of the root to 
be extracted, and whose differences shall also be in arithmetical progression, use 

for expansor —jj-, and then arrange terms with differences themselves differing 

by unity. This may be done by dividing d' , the difference of the expanded fraction 

q' /p' , or q ' —p ' , by n, and making the final interval - „ less than the quotient, 

— after which ascend accordingly. 

Illustrations : 1. Separate f into three factors of above nature : — 

li-=w-6- rhen y- = r 3 ;o' "IT- 4 - 5 ' — =1 - 

Then first interval— 3.5. Therefore : 

17 21.5 27 34 43 54 

13.5' 17' 21.5' ° r 27' 34' 43 

are the desired terms. 

Applying the diagonal process : 

1458 + 1462 + 1462 4382 2191 , 0KQ01Q - „ ^rvuuur 

11 61 + 1161 + 1156 ^WS^IT^ 1 - 25 " 195 - *""'=-0000016. 

2. Separate f into 4 factors : — 
Expansor=16 gives |=|-J ; first interval=6.5. 
[The series of consecutive factors is] 

54.5 62 70.5 _80_ 109_ 124_ l4l_ 160_ 

48 ' 54.5' 62"' 70.5' ° r 96 ' 109 ' 124 ' 141 ' 

3. Separate f into five factors : — 

Expansor-- 12. 5 ; first interval=3, [and the series of consecutive factors is] 
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15.5 19.5 24.5 30. B 37.5 31 39 49 61 75 
12.5' 15.5' 19.5' 24.5' 30.5' ° r 25' 31' 39' 49' 61' 

If the factors are consecutive and in arithmetical progression, then 

All quantities p/q may be represented in n terms, which group in three 
classes as follows : 

V7r) x (t) X \T ) =q/p > where J+ h + l==n - 

ThpT1 n / j (.fc+de)b + k(fa + be)d+l(.ad + bc}f „„ 

1/ </ -j(f c +de)a + Kfa+be)c+l(ad + bc)e v ; ' 

The above is called the three-class process. 

Sometimes a quantity will reasonably resolve into n terms, which group 
in only two classes. Hence the following two-class process : 



£)'*(4)'-« 



, _ j(.ad+bc + 2cd)b + k(ad+bc+2ab)d ._ 

v q ~~~ j(ad+bc+2cd)a + k(ad+bc + 2ab)c ^ '' 



n. / 



Or, the following, simpler but not so good : 

-,/ „ - ic + d)jb+{a+b)kd 
1 q ~~ (c+d)ja+(a+b)kc (0V) - 

Of all these processes the three-class (28) is the most trust-worthy. When 
q does not naturally resolve in such terms, take q/pv which does, and extract 
"l 'v, and multiply results. 
In the consecutive series 

be e b 

— X-j-X X-j-X — =<7, 

a b a e l 

take a new term, of the first expanded q times, that is qb/qa. Then drop the 
term b/a, and call qb=g and qa=f, giving the equation 
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c d f a 

-j-x — x x — x J r-=q- 

be e / 

Now apply the diagonal process, according to the spirit, and not the let- 
ter, of (25), and we have, 

bg+2cf+2de 

2fc/+2ce-M 2 H ~ ] ' q {lbi >- 

Now this result will be found no nearer than the result in (25) but 
the mean of the two will give a close result. It is not, however, a formula of 
value, and I think there are cases where the error of (25) and (251) are both on 
the same side. Hence, their mean would be of no value in particular. 

Process fob Special Roots. 

For Square Root: I Hutton gives — ~— =- — = !-— x — "1 
L. lab \ a & J 

* /„ i (e+d)b + (a+b)d .„ 

2 / 'dbcd+'dabd + b s c+d*a 

./. , (i/^axfe) + (|/-58.d) 

1 H ((/rfXOj+^Xc)' 

or h- b *~B_+ d S~* (33). 

ay/ cd+C]/ ab ' 

If we call r an approximate value of the required square root, then 

^-=,= \pL (34). 

p + r \ p K > 



q 2 +6r 8 q+r 4 
4r(q + r 8 ) ' 



or ^-i-;-^-—*,^* y'q (35) . 



For Cube Roots :- 



7g' + 42q'+80q+2 

2<2S+30q s -|-42q-f-7 ' /q {db >- 



Also, if —x — x -7-=^, 



o c e 
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a / n -)J > d c + b < if+ a 4f 



V ' q ~ l ~acf+ade+bte (37) ' 

This is the best of all cube root processes, and I call it the intenceaving 
process. It has remarkable properties. 

For the | Root : — If ■ — X -j-X — =q ; then 
8 a b c H ' 

!£$-'- V*. or V," (38). 

For the Fourth Root: — Let 

b c d e 
— X -=- X — x -j-=q. 
abed' 

Thpn * > a i (.1>e + cd)b + (ae + bd)e+(ad + bc)d + (ac+b*)e 

VI i (be + cd)a + (ae + bd)b + (.ad+bc)c+(ac+b*)d l >' 

Under-Square Formula. 
U={\/q-\/f)\ U^l/q-Vpy, U a =Wq-Vp)' 

(7) is an under-square formula and, if q is an even square, is precisely 
equivalent to the Jrath root of \/q, by the first correction of formula 1. 
The first correction of (7) is 



n I 



ns+d 5 '~U 5 — U. 

3n bn 

, 2(»»-l) _ « 8 -4 rr ...l«i;. 

ns—d '-U ji — V. 

n on. * 



, , 2(n»-l) rr {E-\)U n(n-2) 
A ls „ "IT , " + i —BT > *-^-*-Z=t- (41 , 

3n 72 w— I 

Still another but complex value of !_3_ is : 

SP 

(2n. s -n) («. 2 -4) J7 4 (n»- 16)(Vg~-» t /p)«x(|/^ + l ) l /p 

ns + a ~ 3n U Gn *12» ' 

d ^-^u («'-4) ^ (w'-iexyg-yp^xd/g+iyp ' ' c )- 

3« 6re *12w 

*It seems that Hit is better in actual practice. 
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For cube root use for coefficients \^and f s instead of V and r 6 s . 

Illustrations : All the equidistant formulae, except these, approach accu- 
racy only when q is near unity. No such restriction binds these [the under- 
square] formulae, especially those which are most developed. And here, let me 
say, is undoubtedly the beginning of a series which I think the Calculi would un- 
fold, and I trust some friend of science will, take the burden of solving it. 
So developed, I am satisfied that all roots of all numbers would be extractable to 
any required degree of accuracy. 

1. Taking (40) extract f/ 40967 

R t— 3 x 4097 + 4095— V- x 3969-yy x 2401 
3 x 4097-4095 - V x 3969- ^ x 2401 

__12291 + 4095- 7056- 667 nearly __ 8663 _, u , 



"12291-4095-7056-667 nearly 

but taking T 5 7 instead of f %— 8^ y-=17 + 



2. Taking *j/4096, we have 

4 x 4097 + 4095-1 x3969-fx 2401 
4 x 4097-4095-f x 3969-* x 2401 = 



3. Take V 4096 - B y (42), 

6 x 4097 + 4095- V- x 3969-$ x 240 1- ^ x *■%«■ X 65 __ 11907.2 

6 X 4097-4095— V + 8969 -f X 2401- ^ X ifx 65 " 3717 ' 

but should=4. 

To Sum a Series : «S=1, r, r 2 , r"- 1 , n terms. 

Let s=ratio + l=r+l, d—r— 1. 

a -— \j:*-w ^- 

s — nd + - ~ 

3s 

If s exceeds 2, this is not accurate enough to be of value. 

Cube Root by Difference Method. 

Take a 3 <q, fc s >q. Call q—a i =A, b 3 -q=B. Then 

a*B + b*A 



aB + hA 



V? (44). 
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Convenient Formula for Roots of 2. 

iOl.+35 + fcl.) 

V2-|- b(r+1 ) (45). 

^ ' 35(r— 11 k ; 

101n-35 + 23!_ii 

6(r+l) 
A Rough Value of q in Terms of E. 

?='= 8=1 < 4b >- 

Logarithms. 

A curious fact, but scarcely a useful one, is : The logarithm of any num- 
ber approaches 

jra "'Vi?x(i/*-i)fl-+(,/2*--i) ? +m ' (47) ' 

wherein K is the logarithm, B the base of the system, m the characteristic of the 
logarithm of q, the quantity. Now if 5=4, the formula becomes 



^-'-^^- + "' (48 >- 



2o 2B ro ») 

If «=«/j>, then J^-H-S-j-JL + m (49). 

If A' 4 =^-, then q-i-B-Cl + g^) (50) 

-Hh^ (51 >- 

Since Napierian logfi=|=Jf AT, and also, since 

44 S'+ 1 

61 2« 

roughly, then 
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44 1 

„ 61 2n . . ,_ n ^ 

l / 9='=44 f roughly (52). 

The term A' also may be understood to mean the logarithm of q to base 4. 
K of 10=f ; then o=4(l + £)=10. 

The error by this method is easily represented by a curve. 

In the fractional processes continually occur the coaddition of fractions, or 
the adding of numerators together, and of denominators also, when the latter are 
not the same. When the fractions, two or more, are "embracing," the results 
are close, and if "harmoniously embracing" then positively accurate geometric 
averages. They are harmonious when the product of the two terms of each frac- 
tion are the same. 

Thus, | and J are embracing, and their co-sum doubled is 2 xjf =2. 727+. 
The real sum is f J=2.767 + . \ 2 , V an ^ V" are embracing, and their co-sum 
tripled is 3 X ||=4.125. The real sum is VoV—4.200. £ and | are harmonious- 
ly embracing. Their co-sum is {£ or I which is their square root and precise 
geometric average. 

Eough addition may be performed by co-addition. Thus, {t + H^If) 
or t \, which doubled differs from |f J as 3817 differs from 3808, or 1 in 423 parts. 

In the foregoing processes for extraction of roots, calculation of logarithms, 
etc., excepting the factor-process, the accuracy increases rapidly as q approaches 
unity. In general they have value only when q lies between .5 and 2. Their 
value should be tested by logarithms. When q exceeds 2, it may generally be 
divided by the nth power of some simple quantity, which will bring it near unity. 
If this be difficult, use the nearest power of 2 (say m) as previously explained. 

Many apparently absurd problems are readily answered by these formulae. 
Thus : Ij the cube root 0/2.5=}! what ™ the 4th root of 3.333 + ? 

Using the 2-class process, formula (29), — 

/19V f i \ > , . J1N (125x3x19) + (645x4) 
(IT) X (tH"* t0 3 « (375xl4) + (645x3) 1 

=\ ieiF==iM?==l-35545, [error here] 



or by (30), 



399 + 132,^ 531 _ 177 
294 + 99 ~~ 393 ~ 131 



=1.35114. 

True answer=l. 35120. 
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Comparison of Fraction. 

It is frequently necessary, or desirable to compare the values of two un- 
wieldly, yet not very unequal fractions, or to ascertain approximately the com- 
parative value of two ordinary fractions quickly, even if only approximately. 
Hence the following : 

Compare b/a with d/c. Suppose them nearly equal to z/y, and divide 
each fraction by it, giving by/az and dy/cz. Subtract unity from each, and we have 

• J a which make equal to , and — - — — which make equal to — -. Then 



az uaz cz vcz 

uazzfvcz (ua=pvc)z 

(uaz) x (vcz) uavcz 

=the relative excess of b/a. over d/c. If u and v are each unity, the process is 
very simple. 

Illustrations: 1. Compare | and T \ with § as measure. 

(l+f)-l=-iV and (l V-Hf)-l=- ft. 

22-16 1 11-8 1 



■=«r. or =i 



22x16 581 ' 11x8x2 581 " 

True answer= s ' e . 
2. Compare || with |-|. Take base=§. 

by _116 dy _ 196 195 + 117 _ 312 __ 1 



az 117' cz ' 195 ' 195x117 22815 73.1+' 



n , 49 + 29 76 1 

By the formula ^^^ _=_. 



1911-1885 1 
True answer=- 



1911 73.5' 

New York, 1S7S—1SS3. 



